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^^ Abstract 

I We introduce a simple agent-based model which allows us to analyze three stylized facts: a fat- 

, tailed size distribution of companies, a 'tent-shaped' growth rate distribution, the scaling relation 

of the growth rate variance with firm size, and the causality between them. This is achieved 
under the simple hypothesis that firms compete for a scarce quantity (either aggregate demand or 
workforce) which is allocated probabilistically. The model allows us to relate size and growth rate 
distributions. We compare the results of our model to simulations with other scaling relationships, 
and to similar models and relate it to existing theory. 
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1 Introduction 



Growth processes have been the object of active research since the ground-laying work of Gibrat 

CN [S] who described growth as a multiplicative stochastic process. By assuming that growth rates are 

^ random variables, and studying the time evolution of the system, he obtained a lognormal distribution 

^^ of company sizes, i.e. a heavy-tailed distribution. The topic has since then received wide attention [4]. 

fK^ Recent empirical studies suggest that the firm size distribution follows a Zipf (power law) distribution 

■^ [TS1[3H]. Various different models exist for the formation of such distributions with power-law tails. 

■ ' In the context of firm growth, the well-known model by Simon [3^ explains a power law for firm size 

f*">) distribution on the basis of a process introduced by Yule [3S] . In this model, the system is constantly 

CO growing, not only in firm size, but also in the number of firms. The exponent of the power law depends 

^"^ on the frequency of new firms, and approaches 1 if this frequency is low. However the need that 

J> the system has to grow continuously for the formation of a power law has led to criticism [3SJ [31] . 

For systems of constant global size, which are the focus of this paper, models exist that explain the 

formation of fat tails by multiplicative stochastic processes. For such systems with multiplicative noise, 

j^ no result of the generality of the central limit theorem exists [37] . However, it has been shown that the 

stationary distribution has a power law tail in the lowest order approximation, when the noise satisfies 

some stringent conditions [TUl HH HU [57] . 

More recently, Stanley et al [19] uncovered two empirical features that an accurate theory on 
firm growth should explain: the growth rate frequencies exhibit an exponential (Laplacian) decay - 
giving rise to a tent shape in logarithmic scale -, and the variance of the growth rate scales with the 
company size n as <j{n) ex n~^ . This means that the simple assumptions by Gibrat and others, who 
assume multiplicative noise to be independent of firm size, are at odds with the data. The empirically 



determined values of /3 (typically ~ 0.2) depend on the studied system. Many papers have provided 
further evidence of these two findings in many growth processes: firm growth [T^ (/? = 0.15), |3H] 
(/3 = 0.18), [30] (/3 = 0.28), [H] (/? = 0.3) (the latter authors also consider bird populations and 
mutual funds), a country's GDP growth [21] (/? = 0.15), citations in scientific journals [3T] (/? = 0.22) 
and the growth rate of crime [32] (/3 = 0.36). The multitude of examples suggest that the process 
generating both a tent-shaped growth rate distribution and a scaling exponent /? 7^ is simple and 
universal. 

A number of sophisticated models giving rise to a tent shaped growth rate distribution have been 
proposed, following different approaches. Botazzi and Secchi [3] predict a tent-shaped (Laplacian) 
growth rate distribution as being the result of a number of abstract shocks drawn from a Polya urn, 
without adressing the question of the variance scaling exponent. In order to obtain the variance 
scaling, many models assume that firms have an internal structure, i.e. they are composed of subunits 
[22l[30l[2^ In Schwarzkopf et al.[2T], the probability that firms' subunits reproduce themselves follows 
a power law. As a result, the aggregated growth rate distribution is also a power law: it is not a 
collective phenomenon but holds already at the individual level. 

Another interesting attempt [21] assumes that the growth rate variance depends on the elements' 
size (which are citations in their case) , and numerically obtains a fat-tailed size distribution. The tent- 
shaped growth rate distribution is however not a result of the model but a hypothesis on individual 
level. 

Most models explaining the tent-shaped growth rate distribution and the variance scaling relation 
do not attempt to explain simultaneously the formation of the fat-tailed firm size distribution. This 
task is tricky: Existing models for power law tails via multiplicative noise assume the growth rate to 
be independent of the firms' size, which seems in conflict with a scaling exponent /3 ^ 0. However, for 
both a power law firm size distribution and a scaling exponent j3 > there is empirical evidence. 

We address this issue in this paper, both theoretically and numerically with a simple agent-based 
model comprising firms and employees. We distiguish between collective phenomena at the firm level 
and at the level of the macroeconomy. In contrast to the models cited above, we investigate the 
hypothesis that the same process accounts for the tent-shaped growth rate distribution, its variance 
scaling relation and the fat-tailed size distribution of firms. 

The paper is organized as follows. In section [2] we introduce the model. In section |3| we discuss 
the size distribution with heavy tails obtained with our model, and give an intuitive explanation in 
terms of a Langevin equation, which is our first main result. In section [4| we present the growth rate 
distribution of the model, which is our second main result, and discuss some consequences of binning 
the data. In section [5| we compare empirically found scaling exponents to our model and point out 
possible extensions of it. In section |6| we discuss a theoretical aspect of the model and compare it to 
further literature, then we conclude and point out extensions. 



2 The model 

Our agent-based model comprises companies and workers, Nc and N^, which are both constant num- 
bers. At each iteration, firms hire workers in order to produce goods, which are then sold. Based on 
their profits, firms adjust the quantity to produce in the next time step. They compete for a limited 
resource, which can be purchasing power of customers, or workforce. Firms demand a discrete quantity 
of this resource, proportional to their size. The proportionality constant reflecting the profits is (l + fi), 
with /i > 0, which in this paper is assumed to be the same for all firmsF] Since /i > 0, firms attempt 
to increase their size, and the stationary state of the system corresponds to full employment. 

A matching algorithm guarantees that the resource is attributed with the same probability to every 
demand: every open position has the same probability of being filled by a worker. This hypothesis, 
which is equivalent to the microcanonical ensemble in statistical physics, and is also used elsewhere in 
the context of growth processes [Ml 122] ■ Interesting dynamics arise if there is shortage of the resource: 
since the iV^, (1 -I- /u) positions are covered at random, the actual number n^^t+i of employees hired by 
a firm i at time i -|- 1, may be smaller than their joboffer n^.t = rii^((l + fi). It can even be smaller than 
the number of employees n^.i in the preceding period, which confronts us with the difficulty that firms 
may receive no worker at all and vanish. The number of active firms would decrease continuously, and 
workers would eventually accumulate in a monopoly, which is avoided in our model by the introduction 
of new firms. To maintain Nc constant, extinct firms are replaced with by new onea^, initialized with a 
number of workers n^f^ drawn from a distribution J-"(n"*^"'). New firms contribute to the total demand 
for workforce in the next period with the quantity (1 -I- /i)n"('". Workers do not stay at their firm but 
are newly placed at every iteration. However our results would hold if only a (constant) fraction of 
workers of each company were newly placed. 

An analytical formulation of our model with a strictly constant number of workers is complicated. 
We relax this constraint, and assume that it is satisfied on average: (^^ Ui^t) = N^, to obtain equations 
that allow to understand the evolution of the system. Then, the probability for an open position to 
be filled can be written as 

where ^^ fii = iV^ + ^- nf^"^ . Since J2i "-"'^"' ^ ^w we may approximate p « (1 -I- /i)~^. Under such 
conditions the probability that a firm with rii workers will receive ki workers in the next period is the 
binomial distribution 

P(A:,|n.)=(^^')p'^'(l-p^-'=^ (2) 

with mean fiip — Ui, which is precisely the number of workers of firm i at the previous time step, and 
variance "fiip^l — p) = ni/i/(l -I- /i)^. The probability that a firm does not get any worker, in which 
case it disappears, is P{0\ni) = (1 — p)"'. 

For large rii, the binomial distributions may be approximated by Gaussian distributions, exhibiting 
the same n-dependence of their variance. Another implementation uses a different rounding method by 
which firms determine their job offer, such that the Gaussian distribution is also a good approximation 
for small firms. It is used in simulations and can be seen as growth of entirely independent subunits. 



^A scenario with heterogeneous fii of firms is pubUshed elsewhere I17| . In that setting, the hinited resource is 
purchasing power of workers, whereas in this paper it is workforce. 

^It is also possible to re-insert firms at a constant rate, in which case a level of active firms will become stationary 
after some time. For this paper, we used a strictly constant number of active firms, since this allows us to present results 
of systems of identical size. 



where subunits are jobs (see appendix Q. The Gaussian approximation of the probabihty for a firm 
of size rii^t to reach size Ui^t+i is written as 

1 



V{n,^t+i = h\n,^t) = 7=e H -, 7 . (3) 

cTjV 27r 



If the growth rate of a firm is defined as 



9t,t = , (4) 



equation pi) yields for the growth rate probabihty density Q (dropping the index t): 



e(ff.|-.) = ^^e'^^'^'-^^\ (5) 

where c = /x/(l + ^)^ for the binomial approximatiorP] The scaling exponent j3 of the growth rate's 
standard deviation, which is defined through 

a{n) ex n~'^ , (6) 

has the value /3 = 0.5. This value for /3 is a general feature of models that explain firm growth as being 
the sum of the growth of independent subunits. In the literature, subunits often represent the sectors 
in which the firm is active [3S1 HU [31]; in our model these are jobs. We further discuss other values 
for /? and the corresponding empirical evidence in section [5] 

3 Analysis of the fat-tailed size distribution 

Simulations of our model lead to a stationary state that exhibits a fat-tailed size distribution of firms 
shown in figure fl] (a), which is however not a power-law. In order to describe it, we briefiy summarize 
the established theory on the formation of power laws, which we illustrate with numerical simulations, 
before describing the results of our model in that context. 

3.1 Formation of power laws by multiplicative noise 

We recall that the evolution of a random variable n is usually described by the following Langevin 
equation 

nt+i ^ 9tnt + ft, (7) 

with gt a random multiplicative noise and ft a random additive noise. Equation ^ has various 
applications in statistical physics, and is used here to determine the probability density fuction p{n) 
in the stationary state. It has been shown in several ways that a multiplicative noise term gt leads to 
the formation of distributions with power law tails for p(n). For certain discrete noise distributions, it 
is even possible to derive p{n) analytically [T1[T2]. For more general noise distributions, p{n) may be 
deduced under some general approximations [101 UHl HH EU • 

In the case of the size distribution of firms p{n), a multilicative noise term gt in equation ([7| means 
that the size of each firm in the system changes by a multiplicative factor drawn from a distribution 



^For the rounding method detailed in the appendix A this constant 



is ^^ 

(1 + M) = 
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Figure 1: Snapshots (taken after 3900 iterations) of (a) the cumulative size distribution and (b) the 
frequency distribution of firm sizes, obtained from a simulation of our model with N^^ = 10^ workers 
and Nc — lO"' companies. The size distribution is fat-tailed but its slope becomes steeper for larger 
firms. However this is not visible in (b) where firms sizes are binned. The slight irregularity at size 
100 in (a) comes from the additive noise: in this setup newly started firms have a size distribution 
J-{n) that is uniform between 1 and 100. 



g, which is independent of the firm's current size. For instance, if g is drawn from a Gaussian for each 
firm, so would be the aggregated growth rate distribution of firms. This is illustrated with a simulation 
shown in figure (pi, where firms' growth rates were drawn from such a Gaussian distribution (instead 
of the growth term from our model). 

The intuitive interpretation for the necessity of the additive term /, already given in section [2J 
holds generally: for a noise distribution where g that can take values < 1, some firms have a nonzero 
probability of shrinking to size 0, i.e. they "die out" . Since the number of companies should also 
be stationary in the stationary state, an additive noise / is needed [lOj . which means for instance to 
re-insert firms having shrunk to size zero, or to add some term to firms whose size is falling below a 
threshold, such as [401139]. which is preventing firms from attaining size 0. 

[10] showed that the size distribution can be derived by expanding its characteristic function in 
the stationary state. At the lowest order approximation, the cumulative distribution P{n) exhibits a 
power law tail: 

P{n' > n) ex n-" (8) 



By re-inserting this expression into the condition for the stationary state, it is possible to derive [T0ll39] 
the condition 

(5"> = 1- (9) 

which can be used to determine the exponent a if the distribution of g is known. We have numerically 
confirmed the exponents a calculated via equation ([9]) for Gaussian G{g), varying both its standard 
deviation a and the mean m around which the Gaussian is centered (see appendix ^ for technical 
details) . 
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Figure 2: Snapshots (taken after 3650 iterations) of the size (cumulative) (a) and growth rate (b) 
distributions, obtained from of a simulation with a (size-independent) Gaussian multiplicative noise, 
of a system with TV^ = 10^ workers and Nc = 10"* companies. The Gaussian distribution appears 
parabolic in logscale. 



3.2 Firm size distribution in our model 

The question now is what changes in this result when the variance of the multiplicative noise term gt 
is size-dependent. Numerical results of our model, which has a size dependent standard deviation of 
Q{g\n), are shown in figured] The stationary state of the firms' size distribution clearly has a fat tail, 
which nonetheless decays faster than a power law. 

We have not succeeded in solving the master equation analytically for our model. However, the 
difference to the case where g is size-independent can be understood through a simple reasoning based 
on relation ([9]): assuming the Gaussian distribution for g, centered around a value m which is slightly 
smaller than 1, Q yields that the smaller the variance of the noise, the larger the exponent a of the 
power law. 

Interpreting a{n) as the slope of the tangent at size n of P{n) in log- log representation, we expect 
that in our model, which has (3 — 1/2, the calculated a{n) increases with n, which is indeed the case, 
as shown in figure [3J For large sizes where additive noise ft does not modify the distribution, this 
explanation seems to hold approximately. It is based on the assumption that since G{g\n) is Gaussian, 
mostly firms of size nt-i close to n will reach size n at time t. Therefore, locally the variance g is 
assumed to be constant such that ^ holds. 

We have performed simulations for several additive noise distributions 7^(71'"^™). The slope of the 
fat-tailed cumulative size distribution depends on the distribution of G{g\n), and it is sensitive to the 
additive term /: the total demand for jobs then varies depending on J2i^?t^- Since we keep the 
number of available workers constant, this leads to a different probability p for a position to be filled, 
and therefore to a different constant c in the growth rate variance. It also modifies the mean m around 
which the growth rate is centered. These factors have been taken into account for the calculation of 
a{n) in figure [s] by equation Q. Obtained values are close to the slopes fitted at certain values of n 
to the cumulative size distribution from the corresponding simulation. 

The result from our model contains a general implication, which is the first of our two main results: 
as soon as a of the growth rate is size-dependent, the fat-tailed size distribution for firms will decay 



faster than power law, provided that the system has constant global size. The size-dependency of the 
growth rate variance might explain why many authors obtain very good fits of the size distribution 
by a lognormal distribution (see |9]|37j[3] and references therein), which also decays faster than power 
law. Simulations of another scaling exponent than in our model are presented in section [5] after our 
second main result in the next section. 
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Figure 3: For Gaussians g, (g"'"') = 1 yields a(n) as a function of firm size n (for scaling exponent 
/3=0.5). Crosses represent the slopes of the tangent of the simulated cumulative size distribution at 
size n. After 3900 iterations, N^ = 10* and N^ = lO'^. 



4 Growth rate probability distribution for specific firm size 
distributions 
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Figure 4: growth rate distribution of the same simulation as figure [T] 

Figure I4] shows that despite the Gaussian Gig\n) in our model, in the stationary state G{g) exhibits 
a tent-shape. In fact, n-dependence makes that the growth rate distribution of G{g\n) is wider for 



small firms and more narrow for big firms. Considering the aggregate Gig), small firms account for 
the fat tails of the tent-shaped Q{g\n), and big firms for the peak. Theoretically, this result can be 
explained as follows: The aggregated growth rate distribution for the Nc firms can be written as 



Si9)^^Y.S{9^\n^), (10) 



or, in the continuous limit: 



/"OO 

g{g)= dng{g\n)p{n), (11) 

Jo 

where p{n) is the firms' size distribution. Since we do not have an analytical expression for the 
simulated size distribution of our model, we evaluate the integral as an approximation for Q(g\n) of 
our model and power-law size distributions of exponent a. For firms' size distributions p{n) ex n~°'~^ 
and scaling exponents /3 = 0.5 the integral yields 

,(,) ^ r nO.^-^^e-^"(-^)^.n ^ ^.2"^ (, - l)--r(l - a) . (12) 



For a = 0, it simplifies to g{g) = r^iYl- ^^ ^^^^ '^^^ ^^^ slopes of G{g) are linear on a double- 
logarithmic scale, i.e. presenting a tent-shape. For a > 0.5, the integral diverges at the origin and can 
only be integrated starting from a finite cutoff value uq, since n~"~^ is not a normalized probability 
density. This includes the case of our model (a ~ 1 and /3 = 0.5). If integrated from a cutoff no. 



equation (12) still yields a tent-shaped g{g), its width depending on the cutoff. Integral (12) can be 
generalized to values of /3 other than 0.5, which is interesting since empirical values are a ~ 1 and 
f3 sa 0.25. In that case, the condition for convergence of the integral becomes a > /3, and for a = 
the expression simplifies to Gig) = Yi3T~^' ^^^ smaller is /3, the less peaked G{g), which is intuitive, 
since if /? = 0, the result is a Gaussian G{g)- 

Notice that the shape of G{g) is not very sensitive to the underlying size distribution: equation 



(12) yields an approximate tent-shaped G{g) even for exponential decay of p{n). This suggests that 



despite the deviation from a power law which our model exhibits in its stationary state, the idea of 



performing integral (11) explains the observed tent-shape well. 

In the literature, the principle of performing this integral has been used in the model by [30) to 
obtain a tent-shaped growth rate distribution of a single firm. Other models [111 UHl 131 [31] do not 
perform the integral and do not clearly distinguish between the growth rate probability at firm level 
and at aggregate level, which they assume to follow a Laplacian distribution. The necessity to perform 



the integral in equation (11) is however independent of the assumed Gig\n). It turns out that the 
functional form of Gig) yields an approximate l/\g — 1| tent-shape for both Laplacian and Gaussian 
Gig\n). If it is integrated from a size cutoff, for low values of /3, Laplacian Gig\n) yield a more peaked 
Gig)- Since the growth rates have typically values close to 1, empirical evidence can often be fitted 
with a Laplacian (centered around 1) and a l/\g— 1|- law equally well. However [30] find that the tails 
of the tent-shape exhibit power law decay rather than exponential decay of a Laplacian, substantiating 
our argument. 

4.1 Artefacts from binning 

As already done for the size distribution, we discuss effects arising from binning data: The empirical 
data in [THl [Ml 1211 [31] show tent-shaped growth rate distributions of different widths depending on 



firms size (or country's size or citations respectively), for which the authors propose a Laplacian fit. 
To do this, firms are grouped according to their size in large logarithmic bins. From the slopes of the 
growth rate distribution on logarithmic scale, cr(n) and its scaling exponent j3 are determined. 

Numerical simulations of our model show that aggregation of growth rates of firms within one 
order of magnitude is sufficient to obtain a growth rate distribution that resembles a tent-shape, in 
spite of the Gaussian Q{g\n) in our model (see figure [5|. The reason is again that p{n) « l/n?. This 
reveals that if an ensemble average is used to determine the shape of Q{g\n)^ its functional form is 
only assessed correctly if the samples from the ensemble have precisely the same size. Already a size 
spectrum of one order of magnitude is enough to exhibit a different form for Q{g\ri). The value of /3 
does not seem to depend on binning. Indeed we find again a ex n^'•^^^ . 
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Figure 5: Histograms for growth rate distributions, (a) Aggregate growth rate distribution for com- 
panies ranging from 1 to 5 • 10^ employees from a system with 10^ employees and a Zipf firm size 
distribution, (b) For Gaussian growth rate probability densities, the growth rate distribution of firms 
within one order of magnitude (its smallest size indicated in the plot) resemble a tent-shape, thus in 
qualitative agreement with the data shown in [12]. (c) The slopes of these approximate l/\g — 1|- 
distributions follow the same scaling relation as the variance of the Gaussian growth rate probability 
densities a{n) ex rT 



-0.5 



The very simple underlying microscopic mechanisms, as well as compatibility with empirical ev- 
idence suggest that Gaussian functions might be a simpler alternative to the commonly assumed 
Laplacian shape for Q{g\n). A further argument is that Gaussian distributions are conjugate priors 
to themselves, so they can be the result of the multiplication Gaussian growth rate probabilities for 



various reasons. 



5 Extension and comparison to empirically found values of (5 

The empirical studies cited in the introduction find smaller values for j3 than our model. These 
are generally explained by firm-intern factors contributing to firm's growth, as by [121 [Ml 1121 HI] ■ 
Intuitively, if the growth of a company was entirely dependant on the decisions of its CEO, there 
would be no reason to assume that a company's size should affect its growth rate variance, so values of 
/3 between and 0.5 are possibly due to a contribution of both internal hierarchical structure as well as 



of a firm's size. An additional factor that can lower /3 are aggregate business fluctuations, which have 
not been considered in the derivation above: If the total available amount of the resource fluctuates, 
so do the means of the Gaussian growth rate distributions which are derived in equation (l5| . This is 
independent of the firm's size. 

As a comparison to these empirical studies, we simulated a system with a Gaussian G{g\n) and a 
scaling exponent /? = 0.25. This /3 is not the result of interactions in the job market as in our model, 
but firm's growth is simulated as Gaussian multiplicative noise where a ex n~°'^^, without specifying 
its microfoundations. 
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Figure 6: (a) cumulative size distribution and (b) growth rate distribution of a simulation with 
Gaussian multiplicative noise with a scaling relation of /3 = 0.25. (after 3000 iterations in a system 
with 10^ workers and 10* companies) 
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Figure 7: (a) Growth rate distribution where companies are clustered into size bins (system with 
10^ companies and 10^ workers), (b) scaling exponent determined from the slopes in (a). It yields 
/3 = 0.26, in agreement with the scaling exponent a oc rj~*'-25 of the Gaussian growth rate probability 
densities with which the system has been simulated. 
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The result of simulations shown in figures l6] and IT] are as expected in between pure multiplicative 
noise /3 = (figure [2| and our model with scaling exponent /3 = 0.5 (figure [I] and |4]). Although the 
cumulative size distribution on log-log-scale appears almost linear in the range n S [100; 1000] (figure 
pi), fitting slopes to it shows the n- variation of a, which is half as pronounced for /3 = 0.25 compared 
to f3 — 0.5. In contrast, the slopes of the tent-shaped growth rate distribution appear less linear than 
for /3 = 0.5. From these simulations and the previous discussions we conclude that there is a tradeoff 
if Gaussian G{g\n) are assumed: the more linear the size distribution is on double-logarithmic scale, 
the less peaked the growth rate distribution is. Whether this Gaussian hypothesis can indeed explain 
data depends on how well a power law fits the cumulative firm size distribution, and how peaked the 
tent-shaped growth rate distribution arising from the same dataset is. The strength of this simple 
hypothesis is that it would allow the explanation of heavy-tailed growth fluctuations as a collective 
phenomenon on aggregate level, without having to assume them on firm level, as e.g. |21| do. Even if 
a Laplacian G{g\n) is assumed, as many authors do, the problem remains that (3^0 will lead to a size- 
dependent approximate power law exponent a{n). Although for Laplacian Gig\n) the n-dependency of 
exponent a is half as pronounced as for Gaussians, it would be visible in data for for size distributions. 

6 Discussion 

Having presented the stationary state of our model (which has /3 — 0.5), as well as simulations of 
systems with /3 = and /3 — 0.25 for comparison, we return to some theoretical aspects of our model. 
Its dynamics can be described on three levels: The noise on elementary (i.e. job) level is the same for 
every element, which can double, vanish, or stay constant (see appendix^]). Iterated more than once, 
this process differs from multiplicative noiserl as soon as a job doubles, in the next iteration both jobs 
follow the dynamics separately. This elementary level allows for the calculation of the evolution of the 
size of the companies, which are the second level. On that level, the growth rate probability density is 
Gaussian, with a size-dependent variance. 



Due to this dependency integral (111 becomes non trivial. The tent-shaped growth rate distribution 
of companies only holds at the level of a whole economy, which is the third level. 

An analogy can be drawn to a physical system, where, due to long-range interactions, the statistics 
of an element and the statistics of the ensemble can diverge. [25j describe this by the term superstatis- 
tics, i.e. statistics of statistics, stating that in physical systems with fluctuations, the Boltzmann factor 
of the system is obtained by integrating the Boltzmann factors of every subsystem over their inverse 
temperatures. The analogy to the model presented here is the following: Instead of a Boltzmann 
factor, the quantity that governs the dynamics of a single firm is Q{g\n), which is depending on n. 
This n-dependence is the result of long-range interactions: the hypothesis that every job is taken with 
the same probability implies that every open position interacts with every available employee. For 



Gibrat's model, performing an integral analogous to equation (11) is not necessary, since the interest 



ing quantity is the growth rate, which is size-independent. For other systems with multiplicative noise. 



an integral analogous to (11) may not be omitted. Physical systems with multiplicative noise, where 
the dynamics depend on the square of a Gaussian variable exhibit Tsallis statistics [25] [14] . Links of 
growth processes to nonextensive thermodynamics are presented in [531 [H] and for the case of the 
tent-shaped growth rate distribution [31jand [M] . 

Finally, we compare our model to two existing models which also exhibit /? — 0.5. fT\ presents a 



*In contrast, the model by |26j use an elementary level which follows a dynamics with multiplicative noise. However, 
since the constituent elements therefore exhibit a power law size distribution, it is more difficult to interpret them. 
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micro-founded model for city growth, in which they study two scenarios, of which one corresponds 
to pure muhiphcative noise (/3 — 0) and yields a power law, and one (which they term linear case) 
in which the growth rate variance has exponent /? = 0.5. For the latter scenario, a size distribution 
other than a power law is derived. This setting differs from our model in that city-dwellers do not 
change the city all at the same time. It corresponds to our model if we simulated it in sequential 
update, a situation where workers drawn at random can change company, and the probability of 
joining a particular company is proportional to its size. In sequential update, Gig\n) would not be a 
distribution (as it is in our model as presented, which corresponds to parallel update) but just one 
value, proportional to the current size. We have simulated these sequential dynamics for comparison. 
If the statistics of G{g\n) are calculated after a given number of changes at the level of workers, very 
similar statistics to our tent-shaped G{g) are obtained. 

The second interesting model with /3 = 0.5 is the widely used model by Simon|55], which has also 
binomial Q{g\n) }39j . Simon's model exhibits a power law for firm size distribution, which is due to 
the fact that the system is constantly growing in employees and in the number of firms. However, it 
has been stated before that the power law is not been found if these two assumptions are not satisfied 

7 Conclusion 

We have introduced a simple agent-based model where the growth rate of a firm is the result of the 
constraints it faces in the markets, in this case, the job market. Its growth rates of firms are size- 
dependent and exhibit the scaling exponent /3 = 0.5, as do other models that describe firm growth as 
the results of independent random processes. In order to keep its size constant, we need to introduce 
new firms to compensate the outgoing flow. The firm size distribution at the stationary state decays 
faster than power law, even though for its probability density in a binned representation it is not 
distinguishable from a power-law. On the basis of this model we show, by using theory on systems 
with multiplicative noise, that for a system of constant size, a scaling exponent /3 7^ implies that 
the size distribution decays faster than power law. The result holds independently of the growth rate 
probability density and for smaller values of /3 7^ in the range of empirical findings. 

The second main result consists of the explanation of a tent-shaped growth rate probability density 
as a collective phenomenon. Our model yields a growth rate probability density for firms that may be 
approximated by a Gaussian. Nevertheless, the aggregate growth rate pobability density of the system, 
for which there is empirical evidence, appears tent-shaped. This tent-shaped form is also found if firms 
are grouped into size bins, and a tent-shaped function is fitted to the growth rates of each bin, without 
the need to assume Laplacian G{g\n). The central idea is to take firm's size distribution into account 
when calculating the growth rate distribution. Small firms account for the 'fat tails', and big firms 
dominate in the peak of the growth rate distribution. For comparison we also show simulations of 
systems of Gaussian multiplicative noise with size-dependent a ex n^'^-'^^. Even for the latter case, 
which is in the range of empirical findings, the growth rate appears tent-shaped. 

A scenario with more economically relevant features has been published by the authors in [17]. 
It may be seen as an application of the results of this paper, where firms compete in addition for 
purchasing power of customers. In that setting, fi is heterogeneous for firms, and can evolve over time. 
The entry and exit process are governed by the financial situation of companies, since they can take 
loans at an interest rate, and at some threshold are declared bankrupt and exit the system. These 
additional features modify the noise term g, but only to an extent that the results presented in this 
paper still hold. 
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Appendices 

A Alternative implementation corresponding to growth of in- 
dependent subunits 

The following setup yields a discrete Gaussian growth rate probability density even for small firms. 
Firms demand on average a quantity of workers fii — n,; (1 + /i), where < /i < 1, which is not 
necessarily an integer number. In our simulations, instead of standard rounding, for every existing 
post j, the joboffer j may be 1 or 2: 

{2 with probability /i 
1 with probability 1 — H 

Then, the joboffer n^ is the sum of the joboffers correponding to the posts of a firm. 

This hi is the offer posted in the job market. Then, the aggregate joboffer J2i '^i is collected. If 
available workforce N^^ is inferior to this offer (which is the case studied here), every open post has a 
probability 

N 1 

P=^^^«T-^ (15) 

L, n^ 1 + M 

of receiving a worker. This attribution on its own would yield a binomial constraint, depending on 
a firm's size, as equation ^ states. However if firms determine their joboffer via equations ( 13 ) and 



(14), the number of received workers follows a symmetric distribution between and 2n, if n was the 
size of the firm in the previous timestep. Combining the probabilities for a job to offer 2 posts with 
the constraints in the job market for each of these joboffers, a single job has a certain probability to 
double, and a certain probability to vanish: 

p{j =2)=q= j-^^ (16) 



(i + m)2 

1 + ^2 


(i + m)2 



Pij-l)-P-j^, (17) 

p(.--0)^,^^ (18) 

(19) 

These probabilities are for single posts. For a firm of size n, the probabilities of receiving k workers 
can be calculated out of these probabilities p and q, in an analogous way as the coefficients of the 
Pascalian triangle are found. It is indeed possible to establish a recursion relation for the coefficients 
C. The probability that a firm of size n will have the size k := 2n — I in the following timestep is given 

by 

Pi2n ~ l\n) = Y. {C{p'-'^-\"-'+'') + Cip'-'^q"-'-^^-')) p'-'= g"-'-^^" (20) 

j=0,n-l+2j>0.l-2j>0 
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In this derivation, the rc-insertion of new firms, which sHghtly increases the variance, has been ne- 
glected. NumericaUy, G{g\n) is less noisy with this rounding methods, compared to the case where 
firms offer precisely (1 + fj-)ni jobs. This rounding method is convenient because it yields a Gaussian 
G{g\n) already for small firms, and therefore more orders of magnitude than for a can be taken into 
account for the confirmation of our theoretical results. 

B Technical details on numerical determination of a 



In simulations, we keep the condition that iV^, remains constant. Whenever a firm goes bankrupt, it is 
replaced by a firm of size n' drawn from distribution !F{n'). The size ^^ n"'^"' of these newly introduced 
firms modifies how many workers are available for existing firms, i.e. the average m around which the 
Gaussian G{g) is centered. If for instance the new firms start with a size drawn from a distribution 
T(n"^'^), with average bigger than the median of already existing firms, less workers are available for 
existing firms, and these will on average shrink, because the Gaussian G{g) is centered around a value 



m slightly smaller than 1. This implies that the larger the average size at restart 



the lower 



the average m, and by integrating equation ([9]), the steeper is the corresponding exponent a. This 
relationship between the additive noise / and the exponent a has been found in for a different additive 
noise term / by [40] and (in continuous time) by [39] : both do not allow firms to shrink below a given 
threshold Smin via an additive term. They find that the bigger Smin, the bigger a. Put differently, 
this means that a declines the bigger the average growth rate is. j39J further states that a declines the 
lower the death rate, and the larger the variance of the growth rate is. We reproduced these results 
both evaluating equation ([9]) as well as in numerical simulations. 
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Figure 8: Snapshots (after 3950 iterations) of the stationary size distribution of a system with N^ — 10^ 
workers and Nc = 10^ firms, with (size-independent) Gaussian multiplicative noise {a = 0.05), centered 
around slightly different values m. The slopes are in agreement with the theoretical values: (a) 
a = -2.15 (b) a = -2.45. 
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